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Abstract 

Crosstalk interference is the limiting factor in transmission over copper lines. Crosstalk cancelation 
techniques show great potential for enabling the next leap in DSL transmission rates. An important issue 
when implementing crosstalk cancelation techniques in hardware is the effect of finite world length on 
performance. In this paper we provide an analysis of the performance of linear zero-forcing precoders, 
used for crosstalk compensation, in the presence of finite word length errors. We quantify analytically the 
trade off between precoder word length and transmission rate degradation. More specifically, we prove 
a simple formula for the transmission rate loss as a function of the number of bits used for precoding, 
the signal to noise ratio, and the standard line parameters. We demonstrate, through simulations on 
real lines, the accuracy of our estimates. Moreover, our results are stable in the presence of channel 
estimation errors. Finally, we show how to use these estimates as a design tool for DSL linear crosstalk 
precoders. For example, we show that for standard VDSL2 precoded systems, 14 bits representation of 
the precoder entries results in capacity loss below 1% for lines over 300m. 
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I. Introduction 

DSL systems are capable of delivering high data rates over copper lines. A major problem 
of DSL technologies is the electromagnetic coupling between the twisted pairs within a binder 
group. Reference [1] and the recent experimental studies in [2], [3] have demonstrated that 
vectoring and crosstalk cancelation allow a significant increase of the data rates of DSL systems. 
In particular, linear preceding has recently drawn considerable attention [4], [5] as a natural 
method for crosstalk precompensation as well as crosstalk cancelation in the receiver. In [2], [3] 
it is shown that optimal cancelation achieves capacity boost ranging from 2x to 4x, and also 
substantially reduces per-loop capacity spread and outage, which are very important metrics from 
an operator's perspective. References [5], [6] advocate the use of a diagonalizing precompensator, 
and demonstrate that, without modification of the Customer Premise Equipment (CPE), one can 
obtain near optimal performance. Recent work in [7], [8] has shown that a low-order truncated 
series approximation of the inverse channel matrix affords significant complexity reduction in the 
computation of the preceding matrix. Implementation complexity (i.e., the actual multiplication 
of the transmitted symbol vector by the preceding matrix) remains high, however, especially 
for multicarrier transmission which requires one matrix-vector multiplication for each tone. 
Current advanced DSL systems use thousands of tones. In these conditions, using minimal word 
length in representing the precoder matrix is important. However, using coarse quantization 
will result in substantial rate loss. The number of quantization bits per matrix coefficient is 
an important parameter that affects the system's performance - complexity trade-off, which we 
focus on in this paper. We provide closed form sharp analytic bounds on the absolute and 
relative transmission rate loss. We show that both absolute and relative transmission loss decay 
exponentially as a function of the number of quantizer bits and provide explicit bounds for the 
loss in each tone. Under analytic channel models as in [9], [10] we provide refined and explicit 
bounds for the transmission loss across the band and compare these to simulation results. This 
explicit relationship between the number of quantizer bits and the transmission rate loss due to 
quantization is a very useful tool in the design of practical systems. 

The structure of the paper is as follows. In section II, we present the signal model for a 
preceded discrete multichannel system and provide a model for the precoder errors we study. 
In section EI, a general formula for the transmission loss of a single user is derived. In section 
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IV we focus on the case of full channel state information where the rate loss of a single user 
results from quantization errors only. Here we prove the main result of the paper, Theorem 14.11 
We provide explicit bounds on the rate loss under an analytic model for the transfer function as 
in [9]. We also study a number of natural design criteria. In section V we provide simulation 
results on measured lines, which support our analysis. Moreover, we show through simulation 
that our results are valid in the presence of measurement errors. The appendices provide full 
details of the mathematical claims used in the main text. 

II. Problem Formulation 

A. Signal model 

In this section we describe the signal model for a precoded discrete multitone (DMT) system. 
We assume that the transmission scheme is Frequency Division Duplexing (FDD), where the 
upstream and the downstream transmissions are performed at separate frequency bands. More- 
over, we assume that all modems are synchronized. Hence, the echo signal is eliminated, as in 
[1], and the received signal model at frequency / is given by 

x(/) = H(/)s(/) + n(/), (1) 

where s(/) is the vectored signal sent by the optical network unit (ONU), H(/) is a p x p matrix 
representing the channels, n(/) is additive Gaussian noise, and x(/) (conceptually) collects the 
signals received by the individual users. The users estimate rows of the channel matrix H(/), 
and the ONU uses this information to send P(/)s(/) instead of s(/). This process is called 
crosstalk pre-compensation. In general such a mechanism yields 

x(/)=H(/)P(/)s(/) + n(/). (2) 

Denote the diagonal of H(/) by D(/) = diag{Yl{f)) and let P(/) = H(/)-iD(/) as suggested 
in [5]. With this we have 

x(/)=D(/)s(/) + n(/), (3) 

showing that the crosstalk is eliminated. Note that with F(/) = H(/) — D(/) we have the 
following formula for the matrix P(/) 

P(/) = (l + D-i(/)F(/))-\ (4) 
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Following [5] we assume that the matrices H(/) are row-wise diagonally dominant, namely that 

ll^iill >> \\hij\\,^i ^ j. (5) 
In fact, motivated in part by Gersgorin's theorem [11] we propose the parameter r(H) 

r(H) = max f ^^f'^^'A (6) 

l<i<N \ \hii\ J 

as a measure for the dominance. In most downstream scenarios the parameter r is indeed much 
smaller than 1. We emphasize that typical downstream VDSL channels are row-wise diagonally 
dominant even in mixed length scenarios as demonstrated in [8]. 

B. A model for precoder errors 

In practical implementations, the entries of the precoding matrix P will be quantized. The 
number of quantizer bits used is dictated by complexity and memory considerations. Indeed, 
relatively coarse quantization of the entries of the precoder P allows significant reduction of the 
time complexity and the amount of memory needed for the precoding process. The key problem 
is to determine the transmission rate loss of an individual user caused by such quantization. 
Another closely related problem is the issue of robustness of linear precoding with respect to 
errors in the estimation of the channel matrix. The mathematical setting for both is that of error 
analysis. Let 

P= (I + D-iF + Ei)-i+E2, (7) 

where 

• El models the relative error in quantizing or measuring the channel matrix H, and 

• E2 models the errors caused by quantizing the precoder P. 

The problem is to determine the capacity of the system, and the capacity of each user, in terms of 
the system parameters and the statistical parameters of the errors. Note that equation ([7]) captures 
three types of errors: errors in the estimation of H, quantization errors in the representation of 
H, and quantization errors in the representation of the precoder P. 

Our focus will be in the study of the effect of quantization errors in the representation of the 
precoder on the capacity of an individual user. Nevertheless, the estimation errors resulting from 
measuring the channel cannot be ignored. We will show that the analysis of quantization errors 
and estimation errors can be dealt separately (see remark 13.21 after lemma 13. 1|) . This allows 
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US to carry analysis under the assumption of perfect channel information. Then, we show in 
simulations that when the estimation errors in channel measurements are reasonably small, our 
analytical bounds remain valid. 

C. System Model 

We now list our assumptions regarding the errors Ei,E2, the power spectral density of the 
users, and the behavior of the channel matrices. 

Perfect CSI: Perfect Channel Information. Namely, 

Ei(/)=0, V/. (8) 

Quant(2~'^): The quantization error of each matrix element of the precoder is at most 2~'^. 
Namely, 

|E2(/),,| <2-^ V/,V2,j. (9) 
DD: The channel matrices are row-wise diagonally dominant. 

r(H(/)) < 1, V/. (10) 

SPSD: The Power Spectral Density (PSD) of all the users of the binder is the same. Namely, 
we assume that for some fixed unspecified function P{f) we have: 

P,(/) = P(/), Vz. (11) 

The main result of the paper. Theorem 14.11 is based on assumptions ([8]), Q, (flOl) . (fTTI) . 

Assumption SPSD can be lifted, as shown in section IXTVl (appendix H). For the sake of clarity 
we present only the simplified result in the body of the paper. 

In order to obtain sharp analytic estimates on the transmission loss in actual DSL scenarios 
we need to incorporate some of the properties of the channel matrices of DSL channels into our 
model. In particular, we will assume 

Werner Channel model: The matrix elements of the channel matrices H(/) behave as in the 
model of [9]. Namely, following [9] we assume the following model for insertion loss 

|H^^(/,£)P = e-2"^^ (12) 
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where ^ is the DSL loop length (in meters), / is the frequency in Hz, and a is a parameter that 
depends on the cable type. Furthermore, crosstalk is modeled as 

|H^^^^(/,^)|2 = mf\R'\f,^)\' (13) 

Here K{£) is a random variable studied in [10]. The finding is that K{t) is a log-normal 
distribution with expectation, denote there ci(^), that increase linearly with I. 

An additional assumption that we will make concerns the behavior of the row dominance of 
the channel matrices H(/, £). 

Sub linear row dominance: 

r(H(/,£)) <7i(£)+72(^)/ (14) 

Where 72(£) = O(V^). 
Remark 2.1: Note that 

|H'^(/,<)| ^ ^" 
The sub-linearity in / follows by studying r(H(£, /)) in terms of random variables behaving 
as K{€}. 

D. Justification of the assumptions 

Perfect CSI is plausible due to the quasi-stationarity of DSL systems (long coherence time), 
which allows us to estimate the channel matrices at high precision. 

Quant(2~'^) is a weak assumption on the type of the quantization process. Informally it is 
equivalent to an assumption on the number of bits used to quantize an entry in the channel matrix. 
In particular, our analysis of the capacity loss will be independent of the specific quantization 
method and our results are valid for any technique that quantizes matrix elements with bounded 
errors. 

Assumption DD reflects the diagonal dominance of DSL channels. While linear precoding 
may result in power fluctuations, the diagonal dominance property of DSL channel matrices 
makes these fluctuations negligible within 3.5dB fluctuation allowed by PSD template (G993.2). 
For example if the row dominance is up to 0.1 the effect of precoding on the transmit powers 
and spectra will be at most IdB. 
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Assumption SPSD (see (fTTI) ) is justified in a system witli ideal full-binder preceding, where 
each user will use the entire PSD mask allowed by regulation. Note that in [3] it is shown that 
DSM3 provides significant capacity gains only when almost all pairs in a binder are coordinated. 
Thus the equal transmit spectra assumption is reasonable in these systems. However we also 
provide in section IXIVI (appendix H) a generalization of the main result to a setting in which 
this assumption is not satisfied. 

Assumption Werner Channel model does not need justification whereas our last assumption, 
sub-linear row dominance was verified on measured lines [3] and can also be deduced analyt- 
ically from Werner's model. In practice, the type of fitting required to obtain 7i(^),72(^) from 
measured data is simple and can be done efficiently. Moreover, the line parameters tabulated 
in standard (e.g., R,L,C,G parameters of the two port model), together with the 99% worst 
case power sum model used in standards [12], provide another way of computing the constants 

III. A General Formula For Transmission Loss 

The purpose of this section is to provide a general formula for the transmission rate loss 
of a single user, resulting from errors in the estimated channel matrix as well as errors in the 
precoder matrix. First, we develop a useful expression for the equivalent channel in the presence 
of errors. This is given in formula (flTI) . Next, a formula for the transmission loss is obtained 
(l30l) . The formula compares the achievable rate of a communication system using an ideal ZF 
precoder as in (HI) versus that of a communication system whose precoder is given by ([7]). This 
formula is the key to the whole paper. Note that we use a gap analysis as in [13], [14]. A useful 
corollary in the form of formula (|34|) is derived. This will be used in the next section to obtain 
bounds on capacity loss due to quantization. 

Let H(/) = D(/) + F(/) be a decomposition of the channel matrix at a given frequency to 
diagonal and non-diagonal terms. Thus D(/) is a diagonal matrix whose diagonal is identical to 
the diagonal of H(/). Also we let SNRi{f) be the signal to noise ratio of the i-th receiver at 
frequency / 

In this formula Pi{f) is the power spectral density (PSD) of the i-th user at frequency /, and 
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nj(/) is the associated noise term. We denote 

<{f) = EMf)\'. (16) 

A. A formula for the equivalent channel in the presence of errors 

We first derive a general formula for the equivalent signal model. The next lemma provides 
a useful reformulation of the signal model in 

Lemma 3.1: The precoded channel Q with precoder as in ^ is given by 

x(/) = D(/)s(/) + D(/) A(/)s(/) + n(/), (17) 

with 

A(/) = (I + D-^(/)F(/))E2(/) - Ei(/)(I + D-^(/)F(/) + Ei(/))"^ (18) 
The proof is deferred to appendix A (section I VIII) . 

Remark 3.2: For our analysis, we will assume that Ei(/) = 0, in which case the formula for 
the matrix A simplifies to 

A(/) = (I + D-^(/)F(/))E2(/). (19) 

The relevance of the formula (fTSi) for the experimental part of the paper (where Ei(/) is not 
assumed to be zero) is explained in the next remark. 

Remark 3.3: In formula (l30l) below we show that the impact of the errors Ei(/) and E2(/) 
on the transmission loss of a user can be computed from the matrix A. Thus, an important 
consequence of the lemma is that the effect on transmission loss due to estimation errors (encoded 
in the matrix Ei(/)) and due to quantization errors (encoded in the matrix E2(/)) can be studied 
separately as they contribute to different terms in the above expression for A. 

B. Transmission Loss of a Single User 

Consider a communication system as defined in (O and denote by B the frequency band of 
the system. We let SNRi{f) be as in (fTSi) and let T be the Shannon Gap comprising modulation 
loss, coding gain and noise margin. Let Ri be the achievable transmission rate of the i-th user 
in the system defined in (O. Recall that in such a system the crosstalk is completely removed 
and therefore 

R,= [ log,{l + r-'SNR,{f))df. (20) 
JfeB 
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Let 

R,{f)=\og,{l + T''SNR,{f)) (21) 

be the transmission rate at frequency / (formally, it is just the density of that rate). Let Ri{f) 
be the transmission rate at frequency / of the i-th user, when the precoder in ^ is used. We 
note that while Ri{f) is a number, the quantity Ri{f) depends on the random variables Ei,E2 
and hence is itself a random variable. Let Ri be the transmission rate of the i-th user for the 
equivalent system in (flTI) . Thus, 



Ri = / R,{f)df. (22) 

JfeB 

By equation (fTTI) . the i-th user receives 

p 

Xiif) = d,,{f)s,{f) + di^i Y. ^^Af>Af) + <f) = ^m(/)(1 + ^^Af)Mf) + N,{f) (23) 

i=i 

where Ni{f) = di^i{f) Yl'j^i^i,jif)-'^jif) + ^iif)- Assuming Gaussian signaling i.e. that all Si{f) 
are Gaussian we conclude that Ni(f) is Gaussian. A similar conclusion is valid in the case of a 
large number of users, due to the Central Limit Theorem. In practice, the Gaussian assumption 
is a good approximation even for a modest number of (e.g., 8) users. Recall also that Gaussian 
signaling is the optimal strategy in the case of exact channel knowledge. Therefore, we can use 
the capacity formula for the Gaussian channel, even under precoder quantization errors. 
Definition 3.1: The transmission loss Li{f) of the i-th user at frequency / is given by 

W) = R,{f)-R,{f). (24) 

The total loss of the i-th user is 

^ Hf)df. (25) 



We are ready to deduce a formula for the rate loss of the i-th user as a result of the non-ideal 
precoder in (flTl) . Our result will be given in terms of the matrix A. Recall that A generally 
depends on both precoder quantization errors E2 and estimation errors Ei. 

Denote by Aj j the {i,j)-th element of the matrix A and let 

Let 
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and 



a.(/) = 5,if)T''SNR,if) = ^\^^Af)\'SNR,if), (27) 

(29) 

r-i^iVi?,(/) + i- ^ 

Note that aj(/) and hence gj(A, /) are independent of the Shannon gap T. The next lemma 
provides a formula for the exact transmission rate loss due to the errors modeled by the matrices 
El and E2. The result is stated in terms of quantities g(A, /) and the effective signal to noise 
ratio, T~^SNR,{f). 

Lemma 3.4: Let H(/) be the channel matrix at frequency / and let Ei,E2 be the estimation 
and quantization errors, respectively as in (|7]). Let Li{f) be the loss in transmission rate of the 
i-th user defined in (|24|). Then 



L,{A, /) = - log2 (1 - A;,(/)(l - g,,(A, /))) , (30) 
where ^^(A, /) is given in (|28] ) and ki{f) is given in (|29l ). 

In particular, if Aj i(/) = —1 the transmission loss is log2(l + r^^S'iVi?j(/)), where SNRi{f) 
is defined in (fT5l) . Finally, if ^i,i{f) 7^ —1 we have 



L,(A, /) < Max (^0, log2 ^ ^^^ ^ J j (31) 
The proof of this lemma is deferred to appendix B (section I Villi) . 
To formulate a useful corollary we introduce the quantities: 

P(f) 

M,{f) = max,^,^ (32) 



tiif) = maxi<j<n\Aij\ (33) 

Corollary 3.5: Let H(/) be the p x p channel matrix at frequency / and let Ei(/), E2(/) be 
the estimation and quantization errors respectively as in ([7]). Let Li{f) be the transmission rate 
loss of the i-th user defined in (j24)) . Assume that tj(/) < 1. Then 

i,(A. /) < log, (i±(P^iMMpM)) (34) 
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Proof: By (|27l) we have 

«^(/) = E §77yI^^-^-^^)I'^^^^^^) ^ M,U)UU)\P - l)^Ari?,(/) (35) 

1 + < 1 + (p - l)M,U)U{ffSNRiU) (36) 
Since |Ai,i(/)| < U{f) we get 

|l + A,,(/)p>(l-t,(/))2 (37) 



Thus by (|28l) we have 



1 _ «.(/) + ! ^ l + (p-l)M.(/)tK/)^iVi?,(/) ^^^^ 



g,,(A,/) |1 + A,,(/)|2- 
Notice that the right hand side is larger than one and using (|3TI) of the previous lemma the proof 
is complete. 

Remark 3.6: We note that under simplifying assumptions, such as assumption SPSD (see ([TT])) 
the above formula reduces to 

L,(A, /) < log2(l + (p - l)t^{f)SNR,{f)) - 2 log2(l - Uif)) (39) 

Under the assumption Perfect CSI, we have A(/) = (I + D"^(/)F(/))E2(/) and since we 
further assumed that the channel matrices H(/) are row-wise diagonally dominant we see that 
A(/) ^ E2(/). Thus, ti{f) ^ 2^'^ and we obtain a bound of the form 

L,(A, /) < log2(l + {p- l)5iVi?,(/)2-2^) - 2 log2(l - 2-^) (40) 
For a statement of a bound of this form see formula (|4TI) of Theorem 14.11 below. 

IV. Transmission rate Loss Resulting from Quantization Errors in the 

Precoder 

In the ZF precoder studied earlier we can assume without loss of generality that the entries 
are of absolute value less than one. Each of these values is now represented using 2d bits {d bits 
for the real part and d bits for the imaginary part, not including the sign bit). We first consider 
an ideal situation in which we have perfect channel estimation. 
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A. Transmission Loss with Perfect Channel Knowledge 

Consider the case where Ei = and the quantization error is given by an arbitrary matrix E2 
with the property that each entry is a complex number with real and imaginary parts bounded in 
absolute value by 2 We will not make any further assumptions about the particular quantization 
method employed and we will provide upper bounds for the capacity loss. We do not assume 
any specific random model for the values of E2 because we are interested in obtaining absolute 
upper bounds on capacity loss. 

The following theorem describes the transmission rate loss resulting from quantization of the 
precoder. 

Main Theorem 4.1: Let H(/) be the channel matrix of p twisted pairs at frequency /, and 
r(/) = r(H(/)) as in ©. Assume Perfect CSI ©, Quaiit(2-'^) ©, SPSD (HU), and that the 
precoder P(/) is quantized using c? > | + log2(l + r{f)) bits. The transmission rate loss of the 
i-th user at frequency / due to quantization is bounded by 

L,id, /) < log2(l + 7(c^, f)SNR,if)) - 2 log2(l - v{f)2~^), (41) 

where 

7(rf,/) = 2(p-l)(l + r(/))22-2^ (42) 

and 

v{f) = V2{l + r{f)). (43) 

Furthermore, suppose d> ^ + log2(l + Vmax) with 

rmax = fnaXfeBir(tl{f)). (44) 

Then the transmission loss in the band B is at most 

/ log2(l + -f{d)SNRi{f))df - 2\B\ log2(l - (1 + w)2-"+°-^), (45) 
JfeB 

where is the total bandwidth, 

-f{d) = 2{l + rmax)\p-l)2-"', (46) 
The proof of the theorem is deferred to section HXl (appendix C). 
We now record some useful corollaries of the theorem illustrating its value. 
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Corollary 4.2: The transmission rate loss Li{A,f), due to quantization of the precoding 
matrix by d bits is bounded by: 

L,(A, /) < log2(l + 7(rf, f)SNR,{f)) - 2 log2(l - v{f)2~^) (47) 

where -f{d, f) = 2{p - 1)(1 + r(/))22-2^ and v{f) = ^2(1 + r(/)). If r(/) < 1, a simplified 
looser bound is given by 

L,(A, /) < 2-'^+3-5 + log2(l + 8{p - l)SNR,if )2~^'') (48) 
For the derivation of the first inequality see (|9T1) in section |Kl The simplified bound is based 

on the estimate — log2(l — z) < 2z valid for < 2; < 0.5. 

The next result is of theoretical value. It describes the asymptotic behavior of Li{d) for very 

large d. 

Corollary 4.3: Under the assumptions of the theorem and assuming that Vmax < 1: 

U{d) = 0(2-^^). 

More precisely, we have 



Remark 4.4: By definition, /(n) = 6{g{n)) if and only if 

We note that for many practical values of the parameters (e.g. SNR{f) = 80dB, d < 20, 
P < 100) the first term in formula (|45l) . involving 2"^^^, is dominant. Since we are interested in 
results that have relevance to existing systems we will develop in the next section, and under some 
further assumptions (e.g. assumptions (fT2l) . (fT3l)). a bound for Li{d) of the form ai2^^'^ + 022^*^ 
where the coefficients 01,02 are expressible using the system parameters. This is proposition 

m 

Ensuring bounded transmission loss in each frequency bin 

We now turn to study the natural design requirement that the transmission loss caused due to 
quantization of precoders should be bounded by a certain fixed quantity, say 0. Ibit/sec/Herz/user, 
on a per-tone basis. Such a design criterion is examined in the next corollary. 

Corollary 4.5: Let t > and let d be an integer with 

d > d{t) (49) 
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d{t) 



With 

' log2(1.25t;(/)2*+Vt/n(2)) if 2* - 1 < f 

0.5\og2{5{p - + rySNRiif) /tln{2)) otherwise 
Then the transmission loss at tone / due to quantization with d bits is at most t bps/Hz. 
Proof: By theorem l4~n the loss at a tone / is bounded by log2 (l + 2~^'^m(/))) — log2((l — 
vif)2-^f). Where «(/) = 2{p - 1)(1 + rfSNR,{f) and vif) = V2{1 + r(/)). 
Using 1 - 2t < (1 - tf we get 

L,{d, /) < log2 (1 + 2-'Mf))) - log2(l - 2v{f)2-'). 
We will show that the inequality 

f i±|4^) < . (50) 



l-2v{f)2-'i 

is satisfied for any d > d(t) as in (|49l ) Let ^ = 2^^^ so that the inequality (l50l) is 



l-2t;(/)^ 

This yields a quadratic inequality of the form 



< 2' (51) 



Az'^ + Bz<T (52) 

with A = u{f), B = 2*+^t'(/) and T = 2* — 1. Using lemma [TOTTI (see section 1X1 - appendix 
D), we see that if > do{t) where 



doit) 



log2(1.25f (/)2*+V(2* - 1) if2*-l< 



t ^ ^ Mi 

AA 



0.5 log2(5(p - 1)(1 + rfSNR,ij)l{2^ - 1)) otherwise 

Then Lj(d, /) < t. But (io(t) < d{t) because 2* - 1 > ln(2)t and the resuk follows. 

Remark 4.6: The qualitative behavior is d{t) ^ ai — log2(t) for very small values of t whereas 
d{t) ^ a2 — 0.51og2(t) for larger values of t. 

B. Applications of the Main Theorem 

We now apply theorem l4~n to analyze the required quantization level for DSM level 3 precoders 
under several design criteria. To that end let Ri be the transmission rate of the i-ih. user (|20l) 
and let Li be the transmission loss of the i-\h user as in (|24)) . The relative transmission loss is 
defined by 

October 9, 2008 DRAFT 



15 



Vi = -^= / W)df/ / Mf)df (53) 

J^i J f£B JfeB 

The design criteria are 

• Absolute/relative transmission loss across the band is bounded. 

• Absolute/relative transmission loss for each tone is bounded. 
Bound on Absolute Transmission Loss 

From now on, we will assume that the transfer function obeys a parametric model as in [9]. 
Thus we assume (fT2l) and (fT3l) . 

To bound the absolute transmission loss we estimate the integral in formula (l45l) of theorem 

sn 

Using the model (fT2l) one can easily see that 

SNR^if) = ^jj-f'"'^ 
Moreover, under the assumption (fT4l) we have a linear bound on the quantity r(H(/, £)) that is, 

r(H(/,£))<7i(^)+72(^)/ 

Where 72) (^) = 0{Vi). Putting these together we can estimate the integral occurring in the 
bound (l45l) and the final conclusion in described in theorem HTSl 

The parameters 7i(^),72(^) enter our bounds through the following quantity. 

« = (1 + -„M)^ + 12(1 + + 240 (J^f) (54) 

Remark 4.7: The quantity pe behaves as 1 + Ci~^^'^ and is close to one for i = 300m. 
We are now ready to formulate one of the main results of this paper: 
Theorem 4.8: Under assumptions Perfect CSI, Quant(2"'^), SPSD, Werner model and sub- 
linear row dominance (see ([8]), (HU), (O, ([HI), dM])) we have 

B 

where 



< ^12-'" + 2-'^+^-^ (55) 



^' = W)^^-^^vi^hp' ^^^^ 

We provide a proof of this result in section IXIlappendix E). 
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Bound on Relative Transmission Loss 

The most natural design criterion is to ensure that the relative capacity loss is below a pre- 
determined threshold. We will keep our assumption that the insertion loss behaves as in the 
model (O, dni). 

Let SNRi = and SNR'- = -S-e~"^ be the Signal to Noise ratios of the i-th user at 



SNR( 

lesi irequencies. we aiso uenoie oy ojv/t = "p"' 'inu oy ojvix j = — jH-, 
Finally, we denote 



the lowest and highest frequencies. We also denote by SNR = and by SNR'i 



c, = ^ \og,iSNR,) + ^ log2(^'.) (57) 

The next proposition shows that Cj provides a lower bound on the spectral efficiency of the i-th 
user. 

Proposition 4.1: Assume that the attenuation transfer characteristic of the channel is given by 
(fT2l) . Then the spectral efficiency is bounded below by 

^Ri > c. (58) 

JD 

The proof is deferred to section |4~n (appendix F). 

Corollary 4.9: Let r]i(d) be the relative transmission rate loss of the i-th user as in (|53] ). 
Assume that the transfer function satisfies (fT2l) and (fT3] ). Then 

Viid) < 02-2d ^ i2-d+3.5 (59) 

where 

^ 4 Pill 

Proof: This is an immediate consequence of the upper bound on the average loss and the 
lower bound on -^Ri. 

Ensuring bounded relative transmission loss in the whole band 

The next corollary yields an upper bound for the number of quantized bits required to ensure 
that the relative loss is below a given threshold. 

Corollary 4.10: Let < r < 1 and let d > d^r) where 



d{r) 



. 5 log2 ( ) otherwise 



Then the relative transmission loss caused by quantization with d bits is at most r. 

The proof is a simple application of the previous bound on the relative transmission loss and 

lemma 110.11 (section |X] - appendix D). 
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V. Simulation Results 

To check the quality of the bounds in theorem 14.11 and its corollaries, we compared the 
bounds with simulation results, based on measured channels. We have used the results of the 
measurement campaign conducted by France Telecom R&D as described in [10]. All experiments 
used the band 0-30 MHz. 

Full band 

For each experiment, we generated 1000 random precoder quantization error matrices E2(/), 
with i.i.d. elements, and independent real and imaginary parts, each uniformly distributed in the 
interval [— 2^'^,2"']. We add the error matrix to the precoder matrix to generate the quantized 
precoder matrix. Repeating this in each frequency we produced a simulation of the quantized 
precoded system and computed the resulting channel capacity of each of the 10 users. Then we 
computed the relative and absolute capacity loss of each of the users. In each bin we picked the 
worst case out of 1000 quantization trials and obtained a quantity we called maximal loss. The 
quantity maximal loss is a random variable depending on the number of bits used to quantize 
the precoder matrices. Each value of this random variable provides a lower bound for the actual 
worst case that can occur when the channel matrices are quantized. We compare this lower bound 
with our upper bounds of theorem 14.11 We have checked our bounds in the following scenario: 
Each user has flat PSD of -60dBm/Hz, the noise has flat PSD of -140dBm/Hz. The Shannon 
Gap is assumed to be lO.Trfi?. As can be seen in figure 1, the bound given by (l45l) is sharp. We 
also checked the more explicit bound (|59| ) which is based on the model (fT2l) . (fT3l) . We validated 
the linear behavior of the row dominance r(H(/)) as a function of the tone / as predicted by 
formula (fl?)) . Next we used (fT2l) to fit the parameter a of the cable via the measured insertion 
losses. The process of fitting is described in detail in [10]. Its value which was used in the bound 
(|59| ) was a = 0.0019. The parameters 71 = 0.1596 and 72 = 3.1729 10^® were estimated from 
the measured channel matrices by simple line fit. The results are depicted in Figure [H 

Single frequency 

The bounds provided for the entire band are results of bounds on each frequency bin. To 
show that our bounds are sharp even without averaging over the frequency band, we studied the 
capacity loss in specific frequency bins. We concentrated on the same scenario as before (i.e. 
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with 10 users), the noise is —lAOdBm/ Hz and the power of the users is —60dBm/ Hz. We 
picked measured matrices H(/i),H(/2), so that SNR{fi) is 40dBm and SNR{f2) is GOdBm. 
As before, we systematically generated an error matrix E2 by choosing its entries to be i.i.d., 
uniformly distributed with maximal absolute value 2^^^+°-^. Next, we computed the transmission 
rate loss using formula (l30l) . By repeating this process = 10000 times and choosing the worst 
event of transmission rate loss, we obtained a lower bound estimate of worst-case transmission 
rate loss. This was compared to the bounds of corollary I4.2[ The results are depicted in figure 
2. Figure 2 uses formula (l47l) . In particular we see that for SNR = 60dB and transmission 
rate loss of one percent, simulation indicates quantization with 13 bits. The analytic formula 
indicates 14 bits. Similarly, when SNR = AOdB, and again allowing the same transmission 
rate loss of one percent, simulation suggests using 10 bits for quantization. The simple analytic 
estimate requires 11 bits. 

The number of quantizer bits needed to assure 99 percent of capacity 

In the next experiment we have studied the number of bits required to obtain a given trans- 
mission loss as a function of the loop length. Figure 3 depicts the number of bits required to 
ensure transmission rate loss below one percent as a function of loop length. We see that 14 
bits are sufficient for loop lengths up to 1200m. Fewer bits are required for longer loops. 

Stability of the results 

In the next experiment we validated that the analytic results proven for perfect CSI are valid 
even when CSI is imperfect as long as channel measurement errors are not the dominating 
cause for capacity loss. To model the measurement errors of the channel matrix H(/), we used 
matrices with Gaussian entries with variance which is proportional to SNR(f). More precisely 
we assumed that the estimation error of the matrix H(/) is a Gaussian with zero mean and with 
variance = tfcttd-ttt, where is the number of samples used to estimate the channel 

H(/) NSNR^(f)' -f 

matrix H(/). For = 1000, we estimated the loss in a frequency bin as the worst case out of 
500 realizations of quantization noise combined with measurement noise. Figure 3 shows that 
as long as the quantization noise is dominant we can safely use our bounds for the transmission 
loss. We comment that the stationarity of DSL channels allows accurate channel estimation. 
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VI. Conclusions 

In this paper we analyzed finite word length effects on the achievable rate of vector DSL 
systems with zero forcing precoding. The results of this paper provide simple analytic expressions 
for the loss due to finite word length. These expressions allow simple optimization of linearly 
precoded DSM level 3 systems. 

We validated our results using measured channels. Moreover, we showed that our bounds can 
be adapted to study the effect of measurement errors on the transmission loss. In practice for 
loop lengths between 300 and 1200 meters, one needs 14 bits to represent the precoder elements 
in order to lose no more than one percent of the capacity. 

VII. Appendix A: Proof of Lemma [3TT] 
In this section we prove lemma [3TT1 

Proof: For simplicity we will omit the explicit dependency of the matrices H(/), D(/), F(/), P(/) 
on the frequency /. We show that 

HP = D + DA, (61) 
with A as above. Indeed H = D(I + D^^F) and thus 

HP = D(I + D-^F)((I + D-^F + Ei)-i + E2). (62) 

Hence, 

HP = D(I + D-^F + El - Ei)(I + D-^F + Ei)"^ + D(I + D-^F)E2. (63) 

Thus, 

HP = D - DEi(I + D^F + Ei)^ + D(I + D F)E2, (64) 
Which proves the lemma. 

VIII. Appendix B: Proof of Lemma [33] 

In this appendix we prove lemma 13.41 

Proof: By equation (flTI) . the i-ih. user receives 

p 

Xi{f) = d,,{f)s,{f) + di^i Y. ^Ms,{f) + n^if) = du{f){l + A,,(/))s,(/) + iV,(/) (65) 
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with Ni(f) = di^i{f) Y^^^i^i,j{f)sj{f) + ni{f). For a large number of users, we may assume 
that Ni{f) is again a Gaussian noise and the transmission rate at frequency / of the system 
described by equation (|65l) will be 

/) - log, ^1 + Ti^^^^p^mi.Af)mMW + \n.(m ) 

Note that this quantity appeared in the main body of the paper just after equation (|22|) where it 
was denoted Ri{f). Dividing both the numerator and denominator by -Pi(/)|(ii,i(/)P we get 

E(Af)- lo. (l + + ^ .... 

R,{A, f) - log, 1 + P,(/) rK(/)P - (67) 

\ z^j^i P^(f)\^t,jyj )\ ^ PdfM,i(f)\^ / 

or 

^.(A./)^.ogil + ia±M4))^) (68, 

\ f'jU J + eSNRiif) J 

where we have defined 

eSNR^if) - - (69) 

and 

'5.(/)=r^Sgl|Ay(/)|2 (70) 



To get the transmission rate loss we denote 

eSNR,{AJ)= (71) 

J + eSNRiif) 

Notice that 

eSNRiif) = eSNRiiO, /) 

By definition (l24l) we have 

L,(A, /) = i?,(/) - R,{A, f) = log2(l + eSNRiif)) - log2(l + eSNR,{A, /)) (72) 
We then have 

r 1 /l + e5iV^^AJ)\ / eSNR,if)-eSNR,iAJ) \ 
L.(A,/) = -log,^ l + e^iVi?.(/) J=-^"g^^^ l + egiVi?.(/) J ^"'^ 

But 

eSNR^if) - eSNR,{A, f) = eSNR,{f) - + ^74^ 

+ eSNR,(f) 
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SO 

.SNR,if) - eSNR4A, /) = ^SNHf W) ^ 1 -1(1 ^ AM))? ,,3, 

+ eSNRiif) 

and finally, 

^^^^ ^^^ eSNR^imif) + 1 - 1(1 + Am(/))P 
e5iVi?.(/) - eSNMA, f) = eSNMf) s.{f)eSNMf) + 1 ^^^^ 

Hence 

L,( A. /) = - log, ^1 - j5^v75(7)Ti — MTrn — j 

where 

a,{f)=6,{f)eSNR,{f) (78) 
and is given in (TTOl) . With the notations (l28l) and (|29l ) we get the formula 

L,(A, /) = - log2 (1 - h{f){l - g,(A, /))) (79) 



To prove the bound we consider two cases. When g(A, /) > 1 we see from equation (1791) 
that Li{A, f) < 0. This clearly indicates transmission gain and the stated inequality is valid. On 
the other hand, if gj(A, /) < 1 we get 

eSNR, .1_^^(^^^))<1_^^(^^^) (80) 



eSNRi + 1 

and using the monotonicity of — log2(l — u) (increasing) in the interval (0, 1), we get 

L,(A, /) < - log2 (1 - (1 - g.(A, /)))) = log2 f — (81) 



g.(A,/) 



and the Lemma is proved. 



IX. Appendix C: Proof of theorem |4JJ 
For the proof of the theorem we need a simple lemma. 

Lemma 9.1: Let A be a complex p x p matrix and define D to be the diagonal matrix with 
Dj j = Aj j for i = 1, Let E be a p x p matrix whose entries are complex numbers with real 
and imaginary parts bounded by 2-"^. Finally, let B = D^^AE. Then \Bij \ < 2-'^+'^/'^{l + r{A)). 

Proof: Let Q = D^^A = I + D"^(A - D). Then we have 

p 

^|Q,,|<l + r(A) (82) 

k=l 
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for all i = 1, ..,p. Therefore 

< 2-^+1/2 J2 l^^^l ^ 2-'^+i/2(l + r) (83) 



p 
fc=l 



fc=l 



Proof of the main theorem 

We first bound the loss Li{f) in a particular tone /. By Lemma [34l we have 



L,(A, /) < Max ( 0, log2 ( ) ) (84) 



where 



,,,A./)^^±Ma^ (85, 

aH/) + 1 

Here A(/) = (I + D(/)-iF(/))E2(/) where H(/) = D(/) + F(/) is the channel matrix at 
frequency / and E2(/) is a matrix whose entries are complex numbers with real and imaginary 
parts bounded by 2^'^. Applying Lemma (19.11) to the matrix H(/) we see that the entries Aij(f) 
are all in a disk of radius v{f)2~'^ around zero. Using r(/) < 5 we obtain = v^(l+r(/)) < 
6^2. Using d > 4 we get 1 - 2''^v{f) > 1 - ^ > 0. 

Thus 

|1 + A,,(/)|2> (l-t;2-'^)2. (86) 
Using the assumption on the PSD of the different users we obtain 



PAf) 

Using Lemma (19.11) we have 



= Y.Wr)\^^M)\"SNMf) = \\Af)\"SNR,{f). (87) 



J2 |A,,(/)r <{P- 1)2-^^+^(1 + r(/))^ (88) 



thus, 



1 + a,(/) < 1 + (p - l)2-''+\l + r{f)fSNR,{f) = 1 + j{d, f)SNR,{f). (89) 
Combining ([85]), (EH) and ((89]) we obtain 

1 ^ l + ^idJ)SNR,if) ^^^^ 



October 9, 2008 



DRAFT 



23 



Note that the right hand side of the above inequality is positive and greater than one. Combining 
IB and (I90l) we obtain 



L,(A, /) < log, ( ^ ^,7^"^' ^^fZ^:^^ ] = log2(l + lid, f)SNR,if)) - 2 log^ll - vif)2-'' 



il-vif)2-^) 

(91) 

Since -fid, f) < 2(1 + W,)'(P-1)2""' and vif) = v^(l+r(/)) < ^2(1 + ^), integrating 
this inequality over f E B we obtain (|45|) and the theorem is proved. 



X. Appendix D: Proofs of Corollary 4.8 and 4.9 

A. A Quadratic Inequality 

In the proof of corollary 4.8 and corollary 4.9 we use the following lemma. 
Lemma 10.1: Let A, B, T be positive real numbers and let 

^(y) ^ I log2(l-255/^) ifT<S 
1 0.51og2(2.5A/T) otherwise 

Then for d > d{T) we have 

^2-2<i ^ 52"'^ < T (92) 
Proof: We let x = 2^'^ and observe that f{x) = Ax^ + Bx is monotone in x > with one 
root of f{x)=T exactly at = ^j^^j^ f^j. ^ > ^^^r^^ ^ \^^^^__^^_) ^e 

have A2-2'^ + 52"'^ = jiX'^) < /(2~'^°) = /(xq) = T. To complete the proof we will show 
that do{T) < d{T). Indeed, 

, , / \ , / 2A( + AAT + B) , 

do T = logo , = logo — — -77^ (93) 

^ ^ ^ \VB^ + 4AT-bJ \ AAT ' 

Thus, 



do{T)=\og,\^{^l + ^ + l)\ (94) 



If we let p = then for p < 1 we have a/1 + p + 1 < 2.5 and this yields the bound 

doiT) < log, (95) 
for T < On the other hand if p > 1 it is easy to see that 1 + ^^1 + p < 2.5y/p thus 
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B A AT \ I A\ 

^o(T)<log2 I— (2.5^-^^)1 =log2 I 2.5^-1 (96) 

Remark 10.2: Note that as T decreases to zero the value of d(T) increases and behaves as 

XL Appendix E: proof of theorem 14.81 
Proof: Using Theorem 14.11 the capacity loss of the z-th user, Li{d), is bounded by 

hid) < [ log2(l + 7(rf, f)^e-'''^)df - 2\B\ log2(l - 2-'^+i-5) (97) 

JfeB 

By assumption, 7((i, /) < 2{p — 1)2"^°'(1 + 71 + 72/)^- To bound the first term we state here a 
simple lemma (for the proof see section IXlIIl - appendix G)). 
Lemma 11.1: Let f{x) = and define 

1 '■^ 



JaAf^) = ^ I log2(l + + hxYf{x))dx (98) 







We have 



^^^^ - ^ + + ) + ^^^^^^ ^^^^ 

We can now finish the proof of the theorem. 

Let a = 1 + 7i, 6 = 72 and /i = 2(p — 1)2^^"', and let J = Ja^ as in the lemma above. From 
(1971) we get 

^Li{d) < J(2(p - 1)2-2'^) - 2 log2(l - 2-^+'-') (100) 

Using the inequality — log2(l — z) < 2z, for z < ^, and the inequality provided by the lemma 
for J(/i) we obtain 

^L^{d) < ^ (2(l + 7iW)^ + 24(l + 7iW)j^ + 240 \og,{l+2{p-l)2-"^f{B))+2~'^'-' 

(101) 

Using log2(l +t) < ln{2)t, the fact that f(B) = -^e""^ and the definition of pe in ^ we 
obtain 
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XII. Appendix F: Proof of proposition 14.11 

Proof: We begin with a bound on the transmission rate of the users. By (|20l) and the model 
(fT2l) we obtain 

i?. = /" log2(l + T-^SNRe~'''^)df > log^ie) [ \n{T-^SNRe-''''^)df (103) 

Jf&B Jf&B 

Thus, 



> B\og,{T-'SNR) - \og,{e) a.^fdf > B\og,(T-^SNR) - -\og,{e)a,BVB (104) 
We notice that this, with SNR' = SNRe''^^'^ implies 

^R, > \og^{V~^SNR)~'^-\og^{e){HSNR)-\n{SNR!)) = ^-\og^{SNR) + ^-\og^{SNR!)-\og^{V), 

(105) 

and the proof is complete. 

Remark 12.1: In practice, the estimation of ai is more reliable than the measurement of the 
transfer function at the edge of the frequency band. Thus, the equivalent form 



^Ri > log,{SNR) - \a,^ 



(106) 



is more reliable. 



XIII. Appendix G: Proof of lemma [UTI] 
In this section we prove lemma lll.li Recall 

1 

J(/i) = — / log^il + ^i{a + hxff{x))dx (107) 
^ Jo 

where f{x) = -^e~°^ 

Lemma: Let a > 1 and 6 > 0. Let M be the maximal value of (a + hxY f{x) in the interval 
[0,5]. We have 

/p^v^ / ah h \ I 

J(/x) < mm I ( + 24— + 240(— )M log^ (1 + /x/(S)) , log^ (1 + M/x) ) (108) 



In particular we have 



2P / , ah f hV\ 
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which is sharp for small values of fi. 

Proof: The inequality J(/i) < log2(l + M/j,) is evident. To get the second bound we compute 
the derivative with respect to // 

1 (a + bxffix) 



J(/x) 



-dx 



Bln{2) ,/o 1 + + bxyf{x) 
Using the lower bound 1 + jj,{a + bx)'^f{x) > 1 + jJ^fix) > 1 + f{B)jj, we obtain 



(110) 



/(/.) < 



a + bxff{x] 



We get 



a2 5Zn(2) (] 

But x"e-^dx = 2 /~ t^^'+^e-^dt = 2(2n + 1)! and hence 



/•B -1 /-oo 

Jo « Jo 

f 

Jo 



f 

Jo 



1 12 
~ Jo 



1 /"°° 240 
x^e-^'^dx < — / x^e-^dx = ^ 
~ Jo 



Thus we get 



P 1 



1 



2a^ 24a6 2406^ 



a2 5/n(2) + 
Integrating this inequality from = to t we obtain 



fj'{n)< 

Jo 



2 P 1 ln(l + t/(5)) 



12ab 1206^ 

+ 



(111) 



(113) 
(114) 
(115) 
(116) 



(117) 



ln{2) al a^B f{B) 
Using the fact that J(0) = 0, we obtain the desired result. 

Remark 13.1: We emphasize that M can be computed analytically. In fact, it is a routine 

exercise to write the maxima M of f{x) in terms of a,b,a. Indeed 

Thus f'{x) = is equivalent to a quadratic equation, and can be solved analytically. Since the 
function f{x) may have at most two critical point, say xi,X2 G [0, oo) we find that 

maxim J{x,)J{xi)J{B)). 
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XIV. Appendix H: Lifting the assumption of equal PSD from the main theorem 

In this appendix we prove a slight generalization of the main result, showing that the assump- 
tion of equal PSD in the binder is not necessary. The resulting bound is similar to that of the 
main theorem 14. 1[ 

To formulate the bound on the transmission loss we introduce the quantities 

PmaAf) = max,{Pi{f)) (118) 

Pminif) = mini;P^(f)^Q{Pi{f)) (119) 

We let p{f) = Pmax{f)/Pi{f). We will say that the PSD satisfies the assumption SPSD(p) 
(or has dynamic range of width p) if we have 

Pmax ( y ) ^ pPmin ( /" ) 

We emphasize that this means that for each / such that Pi{f) 7^ we have 

PmaAf) < PP^U) 

Remark 14.1: In realistic scenarios the number p is limited by the maximal power back-off 
parameter of the modems in the system. 

Theorem 14.2: Assume assumptions Perfect CSI, Quant(2^'^), and SPSD(p). Assume that 
the precoder P(/) is quantized using d > ^ + log2(l + ?"max)) bits. Let H(/) be the channel 
matrix of p twisted pairs at frequency /. Let r(/) = r(H(/)) as in The transmission rate 
loss of the i-ih. user at frequency / due to quantization is bounded by 

U{d, f) < log2(l + j{d, f)SNR,{f)) - 2 log2(l - v{f)2-^), (120) 

where 

7(rf, /) = 2p(/)(p - 1)(1 + r{f)f2-'' (121) 

and 

v{f) = V2{l + r{f)). (122) 

Furthermore, the transmission loss in the band B is at most 

[ log2(l + j{d)SNRi{f))df - 2\B\ \og,{l - (1 + w)2-"+°-^), (123) 
JfeB 
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where is the total bandwidth, and 

-fid) = 2p(l + r^accYip - 1)2-2'^. (124) 
Proof: Only few changes in the proof of theorem 14.11 are needed in order to derive the 
above theorem. In the proof of the main theorem instead of (f87l) we have 

«^(/) = T.%j^\^^M)\'SNR,{f) < J^PimMl'SNR^if). (125) 

The bound on Ajj(/) obtained in (f88l) is valid because our assumptions on the quantization 
are the same as in theorem 14. 1[ Following the same line of reasoning as in equations (|89|)-(|90|) 
yields the bound (|120l) . This, together with the assumption SPSD(p) easily yields (|123l) . 
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Relative Capacity loss at 300m with Perfect CSI 

10^^ F r— . • --r-. 1 I .. I .... I 




No. of quantization bits 

Fig. 1. Relative Capacity loss vs. number of quantizer bits in perfect CSI in a system of 10 users. Integral bound on loss is 
obtained via equation J45b . explicit bound is obtained via ( |59l l and equations J60b . i57\ . 
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Fig. 2. Capacity loss vs. quantizer bits. Perfect CSI in system of 10 users 
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Fig. 3. Number of quantization bits required vs. loop length 
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Relative Capacity loss at 300m: imperfect CSI 
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Fig. 4. Capacity loss vs. quantizer bits. Imperfect CSI in system of 10 users. CSI based on 1000 measurements 
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